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Abstract 

We construct approximate solutions that will describe the last stage of cylindrically symmetric 
gravitational collapse of dust fluid. Just before the spacetime singularity formation, the speed of the 
dust fluid might be almost equal to the speed of light by gravitational acceleration. Therefore the 
analytic solution describing the dynamics of cylindrical null dust might be the crudest approximate 
solution of the last stage of the gravitational collapse. In this paper, we regard this null dust 
solution as a background and perform 'high-speed approximation' to know the gravitational collapse 
of ordinary timelike dust fluid; the 'deviation of the timelike 4-velocity vector field from null' is 
treated as a perturbation. In contrast with the null dust approximation, our approximation scheme 
can describe the generation of gravitational waves in the course of the cylindrically symmetric dust 
collapse. 

PACS numbers: 04.25.Nx,04.30.Db,04.20.Dw 
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I. INTRODUCTION 



The gravitational collapse and spacetime singularity formation are physically very sig- 
nificant phenomena predicted by general relativity. However since extremely high density, 
high pressure and large spacetime curvature will be realized around the spacetime singular- 
ity, all known theories of physics including general relativity could break down there and 
a quantum theory of gravity might be necessary to describe such physical situations. If 
this perspective is right, the appearance of the spacetime singularity predicted by general 
relativity has no rigorous meaning, since general relativity can not describe the gravitational 
processes in the neighborhood of the spacetime singularity. However, it is still important 
to study the spacetime structures around singularities by assuming general relativity, since 
such studies reveal what happens at the entrance of new physics, although there is an issue 
of their observability in connection with the cosmic censorship 

A singularity visible for distant observers is called a globally naked singularity. Naka- 
mura, Shibata and one of the present authors (NSN) have conjectured that the large space- 
time curvature in the neighborhood of non-spherical spacetime singularities can propagate 
away to infinity in the form of gravitational radiation if those singularities are globally 
naked 3|- If this conjecture is true, formation processes of naked singularities are strong 
sources of gravitational radiation and thus might be one of the targets of gravitational wave 
astronomy 

mm 

. However, numerical simulations performed by Shapiro and Teukol- 
sky suggested that the gravitational collapse of a collisionless gas spheroid might form a 
spindle naked singularity in accordance with the hoop conjectured,]^, but little gravita- 
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tional radiation is generated in its formation process[10|]. This seems to be a counter example 
against NSN's conjecture. However it should be noted that the numerical accuracy of the 
simulations performed by Shapiro and Teukolsky seems to be insufficient to give a definite 
statement 12]. 

Motivated by the work of Shapiro and Teukolsky, Echeverria performed numerical simu- 
lations with very high resolution for the cylindrically symmetric gravitational collapse of an 
infinitely thin dust shell [l^]. Although infinitely long cylindrical matter or radiation is un- 
realistic, this system could be a crude approximation of sufficiently thin and long but finite 
length distribution like as the case of spindle gravitational collapse simulated by Shapiro 
and Teukolsky. As Thorne showed, the gravitational collapse of a cylindrical matter with 
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a compact support in the radial direction does not form a black hole horizon (marginal 
surface) Q]. Thus the singularities formed in the simulations by Echeverria must be glob- 
ally naked. Further Echeverria's numerical results showed that the gravitational collapse 
of the cylindrical dust shell causes gravitational wave burst. Thus in this case, the naked 
singularity formation is really a strong source of gravitational radiation. This seems to be 
inconsistent with Shapiro and Teukolsky. 

By the same motivation as that of Echeverria, Chiba investigated the gravitational col- 
lapse of cylindrically distributed dust fluid, which also forms a naked singularity. Then 
Chiba's numerical results imply that the gravitational radiation generated by dust collapse 
is extremely small [l^. This is consistent with the results by Shapiro and Teukolsky and fur- 
ther with Piran's numerical simulations for the gravitational collapse of cylindrical perfect 
fluid with the adiabatic index 4/3 < F < 2; softer equation of state leads to less emission of 

n ' - - 

gravitational waves |13l|. 

In this paper, in order to clarify the reason of the apparent inconsistency between Echev- 
erria's work and the others, and further to understand the generation mechanism of grav- 
itational radiation in the formation of linear singularities, we investigate the dynamics of 
the cylindrically symmetric dust fluid by the analytic procedure. Here we should note that 
it is rather dangerous to accept numerical results, especially these for structures around 
spacetime singularities, without any consistency checks. In this sense, it is also of great 
significance to obtain approximate solutions by completely or almost analytic procedures, 
which could complement the numerical studies. 

A few exact solutions for cylindrical gravitational collapse are known. One is that de- 



scribing the gravitational collapse of null dust 
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which is a cylindrically symmetric 



version of Vaidya solution [17|]. As gravitational collapse proceeds, the speed of collapsing 
matter might approach to the speed of light. Thus the null dust is often regarded as the 
crudest approximation of the last stage of gravitational collapse. However the cylindrical 
null dust does not generate gravitational radiation. Recently, Pereira and Wang constructed 
analytic model of the cylindrical thin shell made of counter rotating dust p articles jisl. This 
exact solution might be useful and physically significant, but unfortunately, it is impossible 
to describe the generation of gravitational radiation near the spacetime singularity [19|. Here 
we would like to comment on Echeverria's analytic work. Echeverria analytically constructed 
an asymptotic solution for the gravitational collapse of an infinitely thin dust shell by using 
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. Although Echeverria's 



the numerical results and further by putting self-similar ansatz 
asymptotic solution might describe the generation of gravitational waves, the behavior of 
the infinitely thin shell might be different from that of finitely distributed masses. Therefore 
in order to see the gravitational wave generation in more general situations analytically, we 
need new exact solutions or new approximation scheme. 

In this paper, we perform a kind of perturbation analysis on the null dust solution by 
assuming that the speed of the collapsing dust is almost equal to the speed of light. We 
stress that the generation of gravitational radiation can be described by this approximation 
scheme. 

This paper is organized as follows. In Sec. II, we present the basic equations of the 
cylindrically symmetric dust system. In Sec. Ill, we derive the equations and solutions of our 
approximation scheme to describe the gravitational collapse with almost the speed of light; 
we explicitly show solutions for thin and thick dust shells, separately. Finally, Sec. IV is 
devoted for summary and discussion of the physical meaning of our solutions by comparing 
those with the previous studies. 

In this paper, we adopt c = 1 unit. 

II. CYLINDRICALLY SYMMETRIC DUST SYSTEM 

In this paper, we assume the whole-cylinder symmetry which leads to the following line 
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ds" = e2(^-'^) {-de + dr^) + e^^dz^ + e-^^R'd^^ (1) 
where 7, ip and R are functions of t and r. Then Einstein equations are 

y = [b!^ - R^y^ I^RR' ^ ^/2^ _ 2RR-tp-ilj' + R'R" - RR' 

-8nG^{R'T\ + RT\)y (2) 

7 = - _ R^y^ ^^RR ^ ^/2^ _ 2RR'^^Ij' + RR" - R'R' 

- SttGv^ {rT\ + R'T\)y (3) 

l-l" = r-^'-^-^V^9T%, (4) 
R-R" = -SttGv^ (t\ + r\) , (5) 



^ + 1^ - ^" - f ^' = (^** + - T\ + T%) , (6) 

where is the determinant of the metric tensor, a dot means the derivative with respect to 
t while a prime means the derivative with respect to r. 

In this paper, we consider dust fluid. The stress-energy tensor is 

T^" = pu^'u", (7) 

where p is the rest mass density and is the 4-velocity field; all these matter variables 
depend on t and r only. By the assumption of the whole- cylinder symmetry, the components 
of the 4-velocity are written as 

u^" = u' (1, -1 + V, 0, 0) . (8) 

By the normalization u^Ufj, — —1, u* is expressed as 

= / (9) 
sJV (2 - V) 

Here we introduce a new density variable D defined by 

D = Re^ 

(2 - V) V{2-V) ^ > 

The components of the stress-energy tensor are then expressed as 

= -e^-^D, (11) 
^g^r^ = e^-^L>(l -V)^ -V^TK, (12) 
y/^T\ = e^-'^D{l - Vf, (13) 

and the other components vanish. The equation of motion V/jT^a = (V/j is the covariant 
derivative) leads to 



d^D = -l(W)' + ^(l-V){2a„(^-7)-V(^-7)}, (14) 
DduV = (1 - V)duD + ^{V{1-V)D}-^ {29„ - ^) - - 7)} , (15) 



2 

2 ^" ' 2 

where u — t — r is the retarded time and du is the partial derivative of u with the advanced 
time V — t + r fixed. 
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To estimate the energy released in the form of gravitational radiation, the C-energy and 
its flux vector proposed by Thornejs^ is useful. The C-energy E is given as 

E:=^{l-e-2'/'(V"/?)(V«/?)}. (16) 

Then the flux vector J'^ associated with E is introduced as 

^Jf^ = ^{drE, dtE, 0, 0). (17) 
Zn 

As shown by Hayward, if the null energy condition is satisfied, ("VaE > for any achronal 
vector and hence < < 1/8 holds outside the causal future of the naked singularity [2^. 
Note that the conservation law V^J^ = is trivially satisfied. Using Einstein equations, 
the energy flux vector J'^ is expressed as 

V^J' = |^|i?i?'(^' + V^'') - 2RRijtlj' - 8nGV^iR'T\ + RT^)^ (18) 

= |^{^^(^' + ^'') - '^RR'i^i^' - SnGy^iR'T't - i?T^)|, (19) 
and the other components vanish. 

III. THE GRAVITATIONAL COLLAPSE OF HIGH SPEED DUST 

Echeverria showed that the speed of a collapsing cylindrical dust shell becomes faster and 
faster and finally almost equal to the speed of light. This seems to be a generic behavior 
near the spacetime singularity even in more general situations, and thus hereafter we focus 
on the situation in which speeds of dust fluid elements are almost equal to the speed of light. 
On the grou nd of physics, such a situation can be well approximated by Morgan's null dust 



solution[lJ]. Then we treat the 'deviation of the 4- velocity of the dust fluid from null' as a 
perturbation and will perform a linear perturbation analysis. 

First, we consider the null dust limit of the timelike dust fluid. The stress-energy tensor 
is rewritten as 

T^u = ^ kfjiy^ (20) 

where 

F = (1,-1 + V, 0, 0). (21) 

In the limit of ^ 0, approaches to the null vector and hence T^^ agrees with that of 
null dust if this limit is taken with D fixed. Here note that the limit of ^ with D 
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fixed leads to vanishing rest mass density p — > 0, i.e., massless field (see Eq. ffTIH) ). Thus 
for < y ^ 1, in order to obtain the approximate solutions for the metric variables, the 
stress-energy tensor could be replaced by that of collapsing null dust. 

Assuming that vanishes initially, the solutions of the complete null dust are easily 
obtained as 

ij = 0, (22) 

7 = Mv), (23) 
R = r, (24) 

SirGDe^^ = (25) 
dv 

where 7R(f) is an arbitrary function of the advanced time This solution was stud- 

n n 

ied by Letelier and Wang|15| and Nolan in detail[16||. If the density variable D does not 
vanish at r = 0, there is the naked singularity; this is not a scalar polynomial singularity 
but intermediate one at which observers suffer infinite tidal force although all the scalar 
polynomials of the Riemann tensor do not diverge 15|. The situation is depicted in Fig.l. 



In this figure, we assume that the density variable D has a compact support with respect to 
the advanced time v, i.e., D > for < f < v^. As mentioned, we will regard this solution 
as a background spacetime. 

We introduce a small parameter e and assume the order of the variables a.sV = 0(e) and 
ip = 0(e). Further we rewrite the variables 7, R and D as 

g7 = 6-^^(1 + 5^), (26) 
R = r{l + 6R), (27) 
D = DBil + 5D), (28) 
and assume that S^, Sji and Sn are also 0(e), where 

We would like to call the perturbative analysis with respect to this small parameter e 'the 
high-speed approximation'. 

The 1st order equations with respect to e are given as follows; the Einstein equations lead 

to 

6^' = SnGe^^DB {6^-ij + 6D- 25„(r5^)} + (r^^)", (30) 



= ^TiGe^^Ds {6^-i) + 6D- 2d,{r6R) -V] + (r^)', (31) 

5^ - 5^" = 0, (32) 

tSr - {r6R)" = IQ-kGc'^DbV, (33) 

- ^" - = ^c'^'^DbV, (34) 

where is the partial derivative of the advanced time v with fixed retarded time u\ the 
conservation law (fT^ leads to 



e 



,7s 



du {5d + 5^ - V) = ^ V— [e-^^DB] + V^'e-^^Ds ^ ; (35) 



the Euler equation ()T5|l reads as 

duV = 0, (36) 

where we have used Eq. fll4|) . Equation ()36p means that V is an arbitrary function of the 
advanced time v, i.e., 

V = V{v). (37) 

This implies that V does not vanish even in the limit to the spacetime singularity. 

Since we solve the basic equations ()3()|) - ()3fi|) as a Cauchy problem, we have to set initial 
data for these equations. Our purpose is to study the generation mechanism of gravitational 
radiation in the formation process of the naked singularity. Therefore we set almost all 
the perturbation variables to vanish at some moment before the formation of the naked 
singularity. More concretely speaking, we set the initial data as 

5^ = Sr = 6r = ip = tp = = at t = ti, (38) 

where is chosen so that the spacetime singularity of the background does not form before 
and at t = ti. Then we may recognize the formation process of the singularity as an origin 
of the perturbation variables generated after t = ti. Here note that these initial conditions 
are consistent with Eq. (jHn|l . 

It should be noted that V and 6^ have to have non-trivial values at t = ti. Since we 
consider the timelike dust, V must not vanish. Due to this constraint and Eq. ()31|) . 6^ has 
non-trivial value at t = ti as 

5^(ti,r) = -S(ti + r), (39) 
where for notational simplicity, we have introduced a new variable S defined by 

S{v):= SirGe'^'DBV. (40) 



For t > ti, the solution of S-y satisfying the initial conditions (|HHjl and is obtained as 

6^{t, r; ti) = I (e(ti - u) r *' +e(n - t^) f *' 1 dxS(t; + x), (41) 

where Q{x) is Heaviside's step function, and we have imposed a boundary condition 6^'\r=o = 
so that the analyticity, or equivalently, the locally Minkowskian nature at the origin r = 
is guaranteed before the formation of the spacetime singularity there. We can also obtain a 
solution of Eq. (jHHj) by the ordinary procedure to obtain the solutions of the wave equation 
in 2-dimensional Minkowski spacetime. Here we should note that tSr has to vanish at r = 
before the singularity formation. The solution satisfying this boundary condition is given as 

2 rt poo POD 

r; ti) = — dr dxSir + x) / — sin(/cr) sin(A;a;) sin {/c(t — r)} . (42) 
7rr Jt; JO Jo k 

We would like to express the solution of Eq. (jH^ in the form 

/ + 00 
dv,S{v,)^{t,r-tuv,), (43) 
-oo 

where by introducing the following function, 

(44, 

\l/ is expressed as 



"^^^'^'^^'^^^ = ^V^^W.) . + (t-.3)/r V T^Ii ■ (45) 

The detailed derivation of the above equation is given in Appendix. Eq. ipKjl can be easily 
integrated and the solution satisfying the initial condition is obtained as 

Sn = ^-S^ + {ti~t) (v^^-2S^ . (46) 

The lowest order of the energy flux vector in the vacuum region is given by 

Hence we need to know the solution for ip to see the energy flux at the future null infinity. 

As mentioned, the symmetric axis r = with non-vanishing Db is the naked singularity 
of the background Morgan spacetime and will also be the naked singularity of the spacetime 



considered here. If S{v) and Vd\n S/dv are functions defined in r > 0, the solutions of 
perturbation quantities PTjl - ljlUj) could be finitely defined in t > ti and thus in the causal 
future of the naked singularity. This means that we could study the structure of the naked 
singularity defined by Eqs. (j4H) - (j46j) . However here note that in general, in order to solve the 
causal future of naked singularities as a Cauchy problem, we need to impose some boundary 
conditions at the naked singularities. In reality, we have imposed the boundary conditions 
for the perturbation equations (jn^ - ()34j) at the naked singularity to obtain the solutions in 
the causal future of the naked singularity; the symmetric axis r = is locally Minkowskian. 

By imposing this boundary condition, we constructed the Green's function of each pertur- 
bation equation and then wrote down the solutions (j4H) - (j43|) . However, after the formation 
of the background naked singularity, it is meaningless to impose the locally Minkowskian 
nature of the axis r = 0, since the axis is singular. Hence we have to be careful to discuss 
the above solutions in the causal future of the naked singularity. For this reason, in this 
paper, we focus on the causal past of the Cauchy horizon only. 

A. Thin Shell 

First, we consider a thin shell, that is, 

S = S,6{v), (49) 

where Sc is a positive constant and 6{v) is Dirac's delta function. We set V = Vc =constant. 
Since the spacetime singularity forms at {t,r) = (0,0), we should choose the initial time 
ti negative. The Cauchy horizon, which is the boundary of the causal future of the naked 
singularity, is given by m = 0. From Eqs.ipUj) and (jlH)), we obtain 

IB = yOiv), (50) 

where we have imposed the boundary condition, 7^ = at r = 0, so that the symmetric 
axis r = is locally Minkowskian before the singularity formation. 

In this thin shell will be shown below, the high-speed approximation breaks down. 

This means that Echeverria's numerical results for the gravitational collapse of a thin dust 
shell can not be recovered by our method. However, the solution of a thin shell is useful in 
constructing solutions for the thick shell case and hence we proceed to study this. 
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From Eq. ()41|) . we obtain the solution for 6^ as 

5,{t, r; to = {0(2ti -u)- Q{v)} . (51) 

From Eq. ip^ . we obtain the solution for 5r in the region of t > ti and of the causal past of 
the Cauchy horizon m < as 

Mt,r;ti) = ^[{t-t-:){Q{v)-Q{u) + Q{u-2U)-Q{v-2t;)} 
2r 

V 7/ 

- -e{u - 2t,)e{v) + -Q{v - 2t,)Q{u)]. (52) 

Substituting Eq.^ into Eq.^, we obtain 

V's(t,r;ti) = 5'c^(t,r;ti,0). (53) 

In order to see the asymptotic behavior in approaching to the Cauchy horizon, we define a 
new variable e as 

t = r{l-e). (54) 

Then the function ^ is expressed as 



r Jo (x + £i)^(a; + £2)(2-£2-2;) 
where we have introduced a new integration variable x = z — f2(t, r; tj, 0) and 



til V 2 

'''' 

The limit of £ ^ 0+ with r fixed corresponds to the limit to approach the Cauchy horizon 
from its chronological past. We can easily see that \E' becomes infinite in the limit of e ^ 0+ 
with r fixed. This divergence of ^ comes from the integral in the neighborhood of x = 0. 
Hence when < e ^ 1, \1/ is dominated by the integral from x = to x = 5x, where 5x is 
some sufficiently small positive constant; 



TT\ r Jo [x + ei)y/x + £2 [V2 



2 \ti\ y^dy 
V r J 

1 [\ti 



TT V r Jo [y^ + £i)vz/ + ^2 



-W^ln£ + 0(1). (58) 
TT V r 
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Thus the asymptotic behavior in the hmit to approach the Cauchy horizon is given as 

S l\t-\ 

Mt,r]U) — > ^W — ln|t-r|. (59) 

n \ r 

Using the above result and from Eq. ()48|) . we find that the energy fiux behaves as 

^a-r (60) 

4:G'K-^{t - r)2 ^ ^ 

The energy /S.E radiated from the neighborhood of the naked singularity is then given by 

= 2n f ^grdt -g^- (61) 



Hence the radiated energy diverges on the Cauchy horizon. This result means that the 
high-speed approximation scheme breaks down there. Using Eq. ()4fij) . we obtain 

5o=^-5, + {U- t) \vc^^^ - 2Se5(t;)| . (62) 

The above equation shows that is infinite and this means that the high-speed approxi- 
mation analysis is not applicable to the thin shell case. 



B. Thick Shell Case 



Let us investigate the gravitational collapse of a thick shell with smooth 7^. Here we 
consider the following S and V\ 

S = Scv^ {v„ - vf e{v)e{v^ - v), (63) 
V = V,v {v^ - v) , (64) 

where v^,, Sc and Vc are positive constants. Also in this case, the Cauchy horizon is located 
at M = 0. Since 6-y and 6r are finite in the causal past of the Cauchy horizon in the thin 
shell case, it is trivial that these are also finite in the thick shell case. Thus hereafter we 
focus on ip and ho- 
using the asymptotic solution (j^Hj) in the thin shell case, we can easily estimate the 
asymptotic behavior of Eq. ()43j) also in the thick shell case. Here we focus on the region 
of sufficiently large r and hence we assume that v^/r is much smaller than unity. In the 
integrand of Eq. (jlK|) . < Vs/r < v^/r <^ 1 holds by the assumption and thus by the same 
procedure to derive Eq. ()58|) . in the limit to the Cauchy horizon e 0+, we obtain 

V'(t,r) ~ _^wlM /" dv^y^^ {v^-v^f\xi(e + —\ — ^-^^^W — Inf— V (65) 
TT V r JO \ r J SOtt V r \ r J 
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where we have also assumed fw/|^i| ^ 1- Thus as expected, ip does not diverge even on the 
Cauchy horizon u = 0. The energy flux is also finite in the same limit; we can easily verify 
that the asymptotic value in the limit to the Cauchy horizon for sufficiently large r is 

^-^^ 5767r3G ^^^^ 

Using Eqs.(|UH|) and Eq. pUj) becomes as 

(5d = - 5^ + 2(ti - t) {VcK - 1v) - S.v'iv^ - vf) 0(t;)0K - v). (67) 

From the above equation, we find that 5z) is finite in the causal past of the Cauchy horizon. 
Together with the result of V', we may say that the solutions of the high-speed approximation 
can describe the asymptotic behavior of the gravitational collapse of the thick shell. 

In the limit of f w with S^vJ' fixed, this model approaches to the infinitely thin shell 
with finite line density. Eqs. ()65|) and (jMj) show that in this limit, and correspondingly 
energy fiux J*" become infinite. This result is consistent with the preceding analysis of 
the thin shell case. The width of the shell corresponds to the time difference of the 
singularity formation; the innermost dust fluid element flrst forms the spacetime singularity 
at t = while the singularity formation time of the outermost fluid element is t = v^. Hence, 
almost simultaneous naked singularity formation will lead to huge emission of gravitational 
radiation. 



IV. SUMMARY AND DISCUSSION 



We analytically constructed approximate solutions for cylindrically symmetric gravita- 
tional collapse of a dust shell with flnite width by the high-speed approximation. This 
approximation scheme could treat the last stage of the gravitational collapse at which the 
collapsing speed is very large due to the gravitational acceleration. We found that thinner 
width of the shell leads to larger amount of gravitational radiation. Hence, if the dust fluid 
of flnite mass collapses into the spacetime singularity almost simultaneously, a gravitational 
wave burst might occur and this is consistent with Echeverria's numerical results. Our re- 
sults also imply that the trajectory of the cylindrical thick shell does not agree with null 
in the limit to the Cauchy horizon and the gravitational emission does not cease there. On 
the other hand, Echeverria's speculation by using self-similar asymptotic solutions is not 
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so in the case of an infinitely tliin shell. Here it should be noted that as shown in Sec. Ill, 
the high-speed approximation breaks down in the infinitely thin shell case and hence our 
results cannot be directly compared with Echeverria's results. Therefore it is a future work 
to reveal the physical reason of this difference between the infinitely thin shell case and our 
present one. 

Our result agrees with Chiba's numerical analysis. In these simulations, the gravitational 
collapse of an infinitesimal portion at r = first forms a massless singularity and then the 
dust fluid in the outer region will accrete on this singularity. In other words, the spacetime 
singularity is gradually formed in Chiba's numerical solutions and hence the gravitational 
waves are not emitted so much in the causal past of the Cauchy horizon. This situation could 
be realized also in the gravitational collapse of the coUisionless gas spheroid in the numerical 
simulations by Shapiro and Teukolskv|lOl|: also in their case, appreciable gravitational waves 
are not generated. 

From our analytic and Echeverria's numerical results, the cylindrical dust collapse could 
be a strong source of gravitational radiation if a finite amount of dust fluid almost simulta- 
neously collapses to form a naked singularity. Of course, the cylindrical distribution is not 
realistic and hence our results do not directly lead to a conclusion that formation processes of 
naked singularities are promising targets of the gravitational wave astronomy. Further stud- 
ies are necessary to obtain more definite implications. However we would like to note that 
recently a possibility of the naked singularity formation process as a candidate of gamma- 
ray burst has been also discussed 2J, |2^. Hence if singularities formed in our universe are 
globally naked, those might have astronomical significance in addition to their importance 
as a laboratory of new physics. 
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APPENDIX A: FORMAL SOLUTION OF tp 



The Green function Gs of the wave equation in 3-dimensional Minkowski spacetime is 
determined by 

(d^ - dl - Ur - ^dl) G,it, r, cp; r, x, = ^{1 - r)5(r - x)6iif - 0), (Al) 

where r and x are radial coordinate [0, +00), and and (p are angular coordinate [0,27r). 
The retarded solution of the above equation is obtained as 

Q(At - Ar) 

G,{t, r, cp; r, x, 0) = ,\ , (A2) 

27rJ(At)2 - (Ar)2 



where 

At := t-T, (A3) 



Ar := yr^ + — 2rxcos(v9 — 0). (A4) 

Using the retarded Green function, the solution for ip is formally written in the form 

ipit^r) := dr dx d(l)S{T + x)Gs{t,r,0;T,x,(l)) 
Jti Jo Jo 

r+00 rt+O roo /•2tt 

dvsS{vs) / dr dx d(f)6{T + x — Vs)Gs{t,r,0;T,x,(f)), (A5) 



ft roo r27r 

r; tj, fs) = dr dx d(j)S{T + x — Vs)Gs{t,r,0;T,x,(j)). (A6) 
Jti Jo Jo 



Hi JO Jo 

Comparing Eq. (|IH|) with Eq. (jA5jl . we find that 

ft roo r27r 

dr dx 

Hi Jo Jo 

Using Eq. ()A2|) . for the region v > Vg and u < Vg, is calculated as 

/•{t—Vs—r) /2+Vs /•+ arccos f2(t,r;T,t)s) 

= TT 

ZTT Jti J — a,rccosQ{t,r;T,Vs) 

1 

X 



y (r — ty — {(r — + 2r(r — Vg) cos0 + r^} 

1 dz r 2{rz+t-vs) dTg 

TT Jn{t,r;ti,v,) y/l - Jti-v, ^(t - ~ - 2{rz + t - Vg)n 



!_ 2\U-Vg\ /•! dz z - n{t,r;tuVs) 

TtV r Jn{t,r;ti,Vs) z + {t — Vs)/r\ 1 — z"^ 

where fl is defined by Eq. f)!^ . In the second equality of Eq. ()A7|) . after introducing new 
integration variables Tg = t — and z = coscp, we change the order of the integrations 
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with respect to Tg and z. Then in the final equahty, the integration with respect to Tg is 
performed. 
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Conical ^ 
Singularity 




FIG. 1: Morgan's cylindrical null dust solution. There is the null dust in the shaded region; 
D{v) > for < V < The dashed hne on i-axis corresponds to the intermediate singularity 
at which an observer suffers infinite tidal force although any scalar polynomials of the Ricmann 
tensor do not vanish. On the other hand, the dotted line on t-axis is the conical singularity. The 
Cauchy horizon is short dashed line of t = r. 
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